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The subject matter of this application is related to the subject matter of co- 
pending U.S. Patent Application Serial No. 09/020,152 filed on February 6, 1998, 
by Konstantine Iourcha, Andrea Nguyen, and Daniel Hung, entitled "FAST 
ADDER/ SUBTRACTOR FOR SIGNED FLOATING POINT NUMBERS", which 
15 is fully incorporated herein by reference. This application claims priority to 
U.S. Provisional Application Serial No. 60/095,995 filed on August 10, 1998, by 
Konstantine Iourcha et al., entitled "RATIO ENGINE FOR GRAPHIC 
INTERPOLATION", which is fully incorporated herein by reference. 

Background of the Invention 

20 1. Field of the Invention 

The present invention relates generally to rendering graphics 
i^Stem, Grid 

^ s yotomo,and more particularly to texture mapping in graphics systems. 
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2. Description of Background Art 

In many graphics applications, three-dimensional images are displayed 
on a computer screen by rendering a large number of graphic primitives, which 
are typically polygons (such as triangles). Primitives are defined using a 

5 number of parameters derived from three-dimensional object models, using 
transformation algorithms, projections, and light calculations. Initially, these 
parameters (or "channel values") are determined for each vertex of the 
primitive. For each vertex, these channel values include: a vertex location (X, 
Y); color values (such as R,G, B) or texture coordinates (U,V), depth of the vertex 

10 location (Z) blending parameter, such as an alpha value, and the like. 

Subsequently, the interior of each primitive is rasterized into a frame buffer 
using texture-mapping or other techniques for assigning color values to points 
within the primitives. Conventional techniques use some form of interpolation 
or incremental calculation to determine on a sequential basis the channel values 

is of the points within the interior of each of the primitives that are rendered. 

Many techniques are known for rendering graphic primitives in such 
graphics systems. For example, Foley et al., Computer Graphics: Principles and 
Practice, 2d. ed., 1990, describes methods for computing edges and for in- 
terpolating values for pixels within the primitives. Conventionally, primitive 

20 edges are determined from vertex points using line interpolation, and scan lines 
between the determined edges are filled with linearly or perspective 
interpolated values for color and other parameters. 

Alternatively, techniques such as "Pixel-Planes" may be used to 
simultaneously determine linear or perspective functions for pixels within the 

25 frame buffer. Parallel linear expression evaluators comprising a binary tree of 
multiply-accumulator nodes are used to evaluate the expressions. See, for 
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example, Fuchs et al., "Pixel-Planes 5: A Heterogeneous Multiprocessor 
Graphics System Using Processor-Enhanced Memories/' in Computer Graphics, 
vol. 23, no. 3, July 1989. This technique has the advantage of improved 
parallelism, but requires expensive dedicated logic for each pixel. 

Another technique, described in Pineda, "A Parallel Algorithm for 
Polygon Rasterization/' in Computer Graphics, vol. 22, no. 4, August 1988, uses 
an edge function. An edge function is a linear function that subdivides a two- 
dimensional plane into three regions: points on the line, points to the "left" of 
the line, and points to the "right" of the line. Boolean combinations of edge 
functions are used to define primitives; for example, three edge functions can be 
combined to define a triangle primitive. 

Once the edges are defined, traversal is employed to fill in the values for 
points within the primitive. Traversal is most often accomplished by 
sequentially scanning rows within the graphic primitive, and using linear 
interpolation to determine a value at each point on each row. See, for example, 
Foley et al., Computer Graphics: Principles and Practice, 2d. ed., 1990, pp. 870- 
71. Also, see Schilling, "A New Simple and Efficient Antialiasing with Subpixel 
Masks", in Computer Graphics, vol. 25, no. 4, July 1991, for antialiasing 
approach employing edge functions to determine distances from relevant edges. 

In conventional scan conversion, processing is divided into two stages, 
namely a setup stage and a rasterization stage. The setup stage receives 
individual primitives that are defined by the relevant values at their vertices. 
The setup stage computes parameters that are needed for performing 
interpolation on an incremental or point-by-point basis. This computation can 
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be an expensive process. In the rasterization stage, primitives are traversed to 
determine values at each individual pixel point. 

Conventional graphics system thus encounter a tradeoff in computational 
load between the setup stage and the rasterization stage. The use of larger 
5 primitives shifts much of the load to the rasterization stage, as fewer primitives 
are used but each requires greater amounts of interpolation. Conversely, the 
use of smaller primitives reduces the amount of interpolation, but necessarily 
involves more computational load in the setup stage due to the increased 
number of primitives. Conventional graphics systems use distinct processor 

20 modules for the two stages, so that computational capacity cannot easily be 
shifted from one stage to the other. Therefore, such systems are often unable to 
efficiently process graphics when the number and size of primitives change. In 
other words, such conventional systems may be optimized for greater 
computational load in one or the other of the stages, but may not easily handle 

25 shifting loads as may occur in dynamic graphic images. 

Another disadvantage of conventional systems is that the rasterization 
stage requires sequential traversal of scan lines within the primitives. In some 
applications, where a value for a given pixel location within a primitive may be 
needed, it may be inefficient to perform lengthy traversal operations of large 

20 portions of the primitive in order to determine the desired value. 

What is needed is a system that combines the setup and rasterization 
stages in such a way that computational load can easily be shifted from one to 
the other without adversely affecting overall performance. What is further 
needed is a system for efficiently interpolating values within graphic primitives 
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without extensive traversal of primitives, and also while avoiding the problems 
and inefficiencies of prior art schemes. 



Summary of the Invention 

In accordance with the present invention, there is provided a system and 

5 method of interpolating values in graphic primitives without requiring separate 

hardware for the setup and rasterization stages. The present invention uses 

direct interpolation to generate a value for any point in a graphic primitive 

without necessarily traversing other portions of the primitive. Unified 

hardware is employed for the graphic rendering operations so that 

w computational load can easily be shifted from one stage to the other. By 

combining the hardware required for the two stages, the present invention 

allows for more flexible resource distribution, and is therefore able to handle 

widely varying categories of images and graphic primitive dimensions. 

In one embodiment in accordance with the invention, a system for 

15 rendering a graphic primitive is provided. The system includes a plurality of 

agents configured to receive information related to the plurality of vertices and 

generate output signals; an arbiter coupled to the plurality of agents and 

configured to receive the output signals and to generate request signals; an 

interpolation engine configured to receive the request signals and generate an 

20 output ratio signal dependent on at least some of the output signals from the 

a 

^ plurality of agents; andja-a-router coupled to the interpolation engine and 

configured to transmit the output ratio signal to an input of at least one of the 
plurality of agents. 
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An interpolation engine is used to interpolate the value of a channel at 
any point on the primitive edge or selected scan line although the practice of the 
present invention is not intended to be limited to a particular interpolation 
engine such as that described below. 

The method of the present invention is generalized so that any desired 
value at a given point can be derived from known values at vertices. Thus the 
known vertices can be described in terms of their channel values, such as 
coordinates (for example, X, Y, U, V, Z, or any combination thereof) or 
parameters (for example, color values), or some combination thereof. Given the 
values of some known parameters at some position, such as a point intersecting 
a line that connects two vertices, the method of the present invention will 
determine, using direct (i.e., non-incremental) perspective or linear 
interpolation, values for unknown parameters at that position. This technique 
provides several advantages. Relatively low setup overhead is incurred, since 
vertices of primitives can be used directly, without additional complex setup 
computations. This results in higher primitive throughput. For example, a 
method of rendering a graphic primitive in accordance with an embodiment of 
the invention includes, deriving a channel value of a first point on a first edge of 
the primitive; deriving a channel value of a second point on a second edge of the 
primitive; and based upon the channel values of the first point and the second 
point, determining a channel value for an interior point located within an 
interior surrounded by the edges of the primitive. 

In addition, this technique facilitates computation on an as-needed basis, 
so that any required computations for other portions of the primitive may be 
deferred until needed. Pixels may be generated in any convenient order, 
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without the sequential limitations of incremental methods. The resulting 
flexibility increases efficiency for implementations of more advanced features or 
techniques. Such methods employ techniques for reordering pixel drawing, and 
therefore are more effective when computation may be performed on an as- 
5 needed or random-access basis. Other visible-surface algorithms can benefit 
from the present invention. See, for example, Foley et al. at pp. 672-95. 

The described technique also avoids several problems associated with 
edge function processing. The present invention avoids needless calculation of 
areas that are not being drawn. Also, since the technique performs interpolation 
10 across horizontal scan lines, sequential horizontal traversal can be made 

extremely efficient. In essence, the present technique facilitates processing for 
points on a random-access basis, but is also capable of providing improved 
efficiency when sequential traversal is used. 

An additional advantage is avoidance of the problem of allocating re- 
15 sources between setup and rasterizing. All of the processing can be handled by 
■=* a single component, so that its utilization and through-put are maximized. 

Whether processing an image with a relatively large number of small primitives, 
or a small number of large primitives, resources are allocated in an efficient 
manner. Shifts in processing load between "setup" and "rendering" are of little 
20 consequence since all functional components may reside in the same piece of 
hardware. 

Finally, the present invention provides smoother gradations of precision 
as coordinate ranges increase. This is useful for processing undipped 
primitives. Since direct floating point computation may be employed, overflow 
25 is rare, and accuracy can be easily controlled. Original vertices can be used for 
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interpolation, even when only a small portion of the primitive (possibly not 
even including the vertices) is being drawn. 

Brief Description of the Drawings 
Fig. 1 is a block diagram of a graphics system for practicing the present 
invention. 

Fig. 2A is a diagram of a primitive to be rendered using the present 
invention. 

Fig. 2B is a diagram of a primitive to be rendered, showing construction 
of line segments for interpolation used in the present invention. 

Fig. 2C is a diagram of a primitive to be rendered, showing additional 
line segments used for interpolation by the present invention. 

Fig. 2D is a diagram of a primitive with a point P to be interpolated 
within the primitive. Fig. 3 is a general block -diagram of an interpolation 
oyotcm according to - thc - prcscnt invention. 

Fig. 2E illustrates a non-degenerative condition for a primitive. 
vFig. 2F illustrates a degenerative condition for a primitive. 

Fig. 4 is a specific block diagram of an interpolation system according to 
the present invention. 

Fig. 5 is a block diagram of an interpolation engine according to the 
present invention. 

Fig. 6 is a schematic diagram showing additional details of a ratio an 
interpolation engine according to the present invention. 
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Detailed Description of the Invention 
Referring now to Fig. 1, there is shown a block diagram of a graphics 
system for practicing the present invention. System 100 is typically imple- 
mented on a computer system such as a personal computer having an Intel 
Pentium™ processor and running the Windows™ 95 operating system from 
Microsoft Corporation. However, those of ordinary skill in the art will realize 
that the practice or implementation of the present invention is not limited to the 
system disclosed in Fig. 1. Rendering engine 101 contains software and 
hardware for rendering primitives into frame buffer 102. Frame buffer 102 is 
typically implemented in an area of random-access memory (RAM) that may be 
optimized for video data storage. Individual memory locations in frame buffer 
102 correspond to defined points on display 103. Data from frame buffer 102 is 
sent to display 103 for output to the user. Display 103 is typically a cathode-ray 
tube, but may be based on some other form of display technology. 

Rendering engine 101 takes input from primitive definitions 104, ren- 
dering parameters 105, and texture maps 106 when rendering primitives. 
Primitive definitions 104 typically include data describing the size, shape, po- 
sition, and other relevant characteristics of graphic primitives. In most systems, 
primitives are polygons such as triangles, so that primitive definitions 104 may 
include coordinates of the polygon vertices, as well as other information, such as 
texture coordinates, for each of these vertices. Primitives may be defined in 
three dimensions using Euclidian coordinates or in four dimensions using 
homogeneous coordinates, and projected onto a two-dimensional plane by a 
known algorithm. In this case, rendering engine 101 or a separate module (not 
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shown) may contain the logic for performing the appropriate transformations 
involved in accomplishing the projection. 

Rendering parameters 105 may include any number of variables and 
parameters that affect the rendering operation. Simulated light sources, 
5 shadows, and other effects may be generated by altering the operation of ren- 
dering engine 101 in a selected manner. 

Texture maps 106 may be used to apply a texture to a primitive being 
rendered by engine 101. Texture maps 106 are typically two-dimensional 
graphic representations of image that is projected onto the primitive to give the 

w illusion that the primitive forms part of a physical surface. 

The present invention typically operates in the context of rendering 
engine 101, though in other embodiments the invention may be applied to other 
components of graphics systems. For illustrative purposes, the description 
given below refers to application of the invention to a graphics system 100 as 

15 shown in Fig. 1. 

As discussed above, primitive definitions 104 typically include the 
channel values associated with the vertices of the primitives. The generalized 
form of the interpolation problem is to determine the channel value of a point 
(such as a pixel) located within a primitive, given the relevant channel values at 

20 the vertices of the primitive. The channel value for a point P within a primitive 
(Fig. 2A) may include but is not limited to any or all of the following channel 
values (or attributes) of the pixel at point P: pixel location (X, Y), color or 
texture coordinates (U, V), and depth (Z) value. One known technique for 
performing this conversion is known as rasterization. As described in Foley et 

25 al. at page 871, various incremental techniques of varying complexity can be 
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used to perform the interpolation. These techniques typically involve scanning 
the primitive on a line-by-line basis, and applying incremental addition to each 
successive pixel in order to derive a pixel value. Considerable processing time 
and setup are required to perform these operations, and random access to any 
given point in the primitive is not easily accomplished without traversing the 
primitive up to that point. 

Referring now to Figs. 2A through 2C, there is shown a series of con- 
ceptual diagrams illustrating the general principles of the present invention. 
Fig. 2 A shows a portion 200 of frame buffer 102 containing a primitive 201. For 
illustrative purposes, only one primitive 201 is shown. Primitive 201 is a 
triangle, though any shape may be used. Primitive 201 is defined by its three 
vertices 202, 203, 204, labeled Vi, V2, and V3, respectively. Accordingly, three 

sides 206, 207, 208 are defined, which form the edges of primitive 201. Point 
205, designated as P, is selected as the point (or pixel) to be rendered. 

In order to determine a channel (or pixel) value for point P from the 
relevant channel values at vertices Vi, V2, and V3, the present invention 

performs direct interpolation as follows. Referring now to Fig. 2B, there is 
shown the same primitive 201 in frame buffer portion 200. A line segment 209 
has been added, defined as a horizontal line segment that intersects point P. In 
alternative embodiments, the line segment may have an orientation other than 
horizontal; however the horizontal line segment as shown is advantageous since 
it facilitates sequential traversal along scan lines. Additional points 210 and 211, 
designated as Pi and P2, respectively, are defined as the points of intersection of 

line segment 209 with any two sides of primitive 201. 
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Referring now to Fig. 2C, there is shown the same primitive 201 in frame 
buffer portion 200. Various distances between points are indicated and labeled, 
as follows: 



» A: 


distance between Vi and Pi; 


> B: 


distance between V2 and Pi; 


» C: 


distance between Vi and P2; 


> D: 


distance between V3 and P2; 


> E: 


distance between P and Pi; and 


> F: 


distance between P and P2. 



Derivation of a pixel value for point P is done by three interpolation 
operations. First, a value for point Pi is derived by interpolation between the 
values for points Vi and V2, as shown in Equation 1 below. Point Pi is a point 
that intersects a primitive edge that is bounded by points Vi and V2. The 
interpolation operation (for linear interpolation) is described by the formula: 

channel value(Pi) = 

(B*channel value of (Vi) + A* channel value of (V2)) / (A+B) 

(Eq. 1) 

Equation 1 describes an interpolation approach that defines a channel 
value of a point along an edge of a primitive as a ratio of the point's distance 
from each of the vertices that define the ends of the primitive edge. Thus, if the 
channel value sought for Pi is its coordinate location, the ratio of Pi's distance 
from Vi and V2, respectively, is calculated and then used as Pi's coordinate 
location. As described more fully below, in one embodiment an interpolation 
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system is employed to perform this calculation. However, it is understood from 
the teachings of the present invention that these calculations may be performed 
other suitable systems. 

Similarly, a channel value for point P2 is derived by interpolation 
5 between the values for points Vi and V3, as follows: 

channel valued) = 

D*channel value(Vi) + Cchannel value(V 3 )) / (C+D) 
(Eq- 2) 

Once values for Pi and P2 have been determined, a third interpolation 
10 between the two is performed in order to derive a value for point P: 

channel value(P) = (F*channel value(Pi) + E*channel value^)) / (E+F) 
(Eq. 3) 

It can be seen that for successive points on the same horizontal span, only 
the third interpolation need be performed. Thus, the present invention offers 
15 added efficiency when sequential traversal along horizontal scan lines is used. 
Although this is advantageous, it is not a required feature of the present 
invention, and in alternative embodiments, points Pi and P2 need not lie on the 

same horizontal line segment. 

Reference is now made to Fig. 2D to further show an example of 
20 performing direct linear interpolation in accordance with an embodiment of the 
present invention. Each of the vertices Vi, V2, and V3 will have channel values 
in the coordinates (X, Y, U, V, Z). The channel values for each vertex is as 
follows: vertex Vi = (X, Y, U, V,Z) = (Vix, Viy, Viu, Viv), vertex V 2 = (V2X, V 2 y, 
V2U, V2v), vertex V3 = (V3X, V3Y, V3u, V3V). 



25 
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(Each of the vertices may also include a channel value Z although this is 
not intended to be limiting in any way. The channel value Z is not required 
when performing the method of direct linear interpolation but is required when 
performing the method of direct perspective interpolation. As shown further 
below, the method of direct perspective interpolation uses the channel value Z 
as a perspective correction coordinate.) 

The Y coordinate for a particular point Pi will be determined by the API 
or equivalent programming interface. For example, point Pi may have a Y 
coordinate of Piy = 2. The channel values of point Pi are as follows in this 
example. Point Pi = (X, Y, U, V) = (Pix, Piy, Piu, Piv). A linear ratio Rpiy may be 
calculated based on the following: Rpiy = (Viy - Piy)/(Viy - Viy). The X 
coordinate Pix of point Pi is then calculated based on the following: Pix = 
Rpiy*V2x + (1 - Rpiy)*Vix. The U coordinate Piu of point Pi may be calculated 
based on the following. Piu = Rpiy*V2u + (1 - Rpiy)* Viu. The V coordinate Piv of 
point Pi may be calculated based on the following. Piv = Rpiy*V2v + (1 - Rpiy)* 
Viv. 

The channel values of point P2 are as follows in the above example. Point 
P2 = (X, Y, U, V) = (P2X/ P2Y, P2U/ P2v). A linear ratio Rp2Y may be calculated 
based on the following: Rp2y = (Viy - P2Y)/ (Viy - V3Y). The X coordinate P2X of 
point P2 is then calculated based on the following: P2X = Rp2y*V3x + (1 - 
Rp2y)*Vix. The U coordinate P2U of point P2 may be calculated based on the 
following. P2U = Rp2y*V3u + (1 - Rp2y)* Viu. The V coordinate P2V of point P2 
may be calculated based on the following. P2V = Rp2y*V3v + (1 - Rpiy)* Viv. 

The channel coordinates of point P is represented as (X, Y, U, V) = (Px, Py, 
Pu, Pv). In the example above, the Y coordinate of point P is Py = 2. The U,V 
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coordinates of point P can be then determined as follows. A linear ratio Rpx is 
determined based on the following: Rpx = (P2X - Px)/ (P2X - Pix). The X 
coordinate Px of point P is calculated as follows: Px = Rpx*P2x + (1 - Rpx)*(Pix). 
In the alternative, X coordinate Px may be provided by the API or equivalent 

5 programming interface. The U,V coordinates Pu,Pv of point P are calculated as 
follows: Pu = Rpx*Piu + (1 - Rpx)*(P2u), and Pv = Rpx*Piv + (1 - Rpx)*(P2v). 

In an alternative embodiment, the API may provide a constant X 
coordinate value. For example, for point Pi, Pix = 2. As a result, a ratio Rix may 
be calculated based on the following: Rix = (Vix - Pix)/ (Vix - V2X). The Y,U,V 

10 coordinates of point Pi are calculated based on the following equations: Piy ~ 
Rix*V 2 y + (1 - Rix)*Viy, Piu = Rix*V 2 u + (1 - Rix)*Viu, and Piv = Rix*V 2V + (1 - 
Rix)*Viv. Similarly, the coordinates P 2 y, P 2 u, and P 2 v for point P 2 may be 
calculated based on the approach similarly described above. The coordinates 
Py, Pu and Pv of point P are then calculated based on an approach as similarly 

25 described above. 

In an alternative embodiment, perspective interpolation is performed in 
lieu of linear interpolation. To perform perspective interpolation, the system 
considers primitive 201 to have a three-dimensional position, so that a depth 
coordinate is available in addition to the conventional two-dimensional X and Y 

20 coordinates. Distances A through F are weighted by a perspective correction 
parameter. In accordance with a preferred embodiment of the present 
invention, the perspective correction parameter used is the channel value Z 
although this approach is not intended to be limiting in any way. The inverse 
form of the depth correction variable W may also be used. In one embodiment, 
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the system of the present invention can be set to one of two modes, for linear 
and perspective interpolation respectively. 

In one embodiment, two pipes operate in parallel: a ratio pipe and 
blending pipe. In an alternative embodiment, a third pipe 266 may be 
5 optionally implemented as shown in Fig. 6. The ratio pipe takes as its inputs 
values for the following variables which are also shown in Fig. 6: 

• X: value of an independent variable, such as physical position, at the 
point being processed; 

• Xo, Xi: endpoint (vertex) values for the independent variable; 
20 • Zo/ Zv endpoint (vertex) values for depth (perpendicular to the 

projection plane); 

• Uo, Ui: endpoint (vertex) values for a dependent variable, such as 
color or texture coordinates, to be interpolated; 

• Vo, Vi: endpoint (vertex) values for a dependent variable, such as 
is color or texture coordinates, to be interpolated; 

In the implementation described below, parameter blenders for color and 
alpha channels (R, G, B, A) are implemented outside the engine in a pixel 
processor (not shown) in a conventional manner. Input Z-values are assumed to 
be positive, and the interpolation is always done within the line segment. 
20 Degenerative line segments (where both end points are the same), are 
permissible. 

In one embodiment, main pipe input is provided in the form of the 
following signals, which may be in an integer or floating point format: 



Variable 



Description 



Format 
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R_X_0 


input variable 0, might be any of X, Y, U, V, 
Z 


1.8.15=24, float 


R_X_1 


input variable 1, might be any of X, Y, U, V, 
Z 


1.8.15=24, float 


R_Z_0 


input (weight) variable 0, usually Z 


1.8.15=24, float 


R_Z_1 


input (weight) variable 1, usually Z 


1.8.15=24, float 


R_X 


requested value of the variable to interpolate 
to 


1.8.15=24, float 


R_U_0 


input variable 0, might be any of X, Y, U, V, 
Z 


1.8.15=24, float 


R_U_1 


Input variable 1, might be any of X, Y, U, V, 
Z 


1.8.15=24, float 


R_V_0 


input variable 0, might be any of X, Y, U, V, 
Z 


1.8.15=24, float 


R_V_1 


input variable 0, might be any of X, Y, U, V, 
Z 


1.8.15=24, float 



In one embodiment, the interpolation engine generates the following 
signals as output: 



5 


Interpol 


ation Engine Output 




Variable 


Description 


Format 




R_U 


Output variable, might be any of 


1.815=24 or float 






X,Y,U,V,Z 
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R_V 


Output variable, might be any of 
X,Y,U ( V,Z 


1.8.15=24 or float 








R_R 


Output ratio 


8 unsigned 



The interpolation engine performs computations as follows. Note that 
the prefix "R_" is omitted in the variables in the calculations below. 



Interpolation calculations are performed as follows: 

IF (Xi - X 0 != 0) { 

U=(Xi-X) ((U 0 -Ui)/(Xi-X 0 >) + Oi 
V=(Xi-X) ((V 0 -Vi)/(Xi-X 0 )) + Vi 

} 

ELSE { 

U = Ui 
V = Vi 

} 

It is noted, however, that if Xi - Xo = 0, other values may also be set for U and V. 
The use of IF/ THEN conditions in the interpolation calculations immediately 
above describes one embodiment for contending with variable instability, such 
as when a divide-by-zero conditions occurs when computing a ratio, and is not 
intended to limit the present invention in any way. A divide-by-zero condition 
is commonly referred to as a "DEGEN" or degenerate case. 

Fig. 2E illustrates a non-DEGEN or non-degenerative case. For example, 
in the non-DEGEN case, the X coordinates of vertex Vi and Vo are not equal (i.e., 
Vox - Vix does not equal zero (0)). The U,V channel values are then calculated 
using the approaches mentioned above. 
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Fig. 2F illustrates a DEGEN case in which Vox = Vix. Hence, when 
computing the ratio values described above, a divide-by-zero condition is 
present. Thus, in the "else" command in the interpolation calculations 
immediately above sets U = Ul (i.e., Piu = Vix). As another alternative, the 
"else" command can set the following condition in a DEGEN case: Piu = Vox . 

Calculations for the perspective ratio for use in perspective interpolation 
are as follows: 

IF { ( ( (Xi - X)Z 0 ) - ( (X 0 - X)Zi) ) !=0) 

PR= ( (Xi-X) Z 0 ) / ( ( (Xi-X) Z 0 ) ~ ( (X 0 -X) Zi) ) 

ELSE 

PR = 0 
where Z 0/ Z,!=l. 

The perspective ratio, PR, may then by used in lieu of a linear ratio to interpolate 
a corresponding channel value. For example, referring again to 2D, the U coordinate for 
point P is interpolated using the following equation: P u = PR (P lu ) + (1-PR) (P 2u ). 
Similarly, the V coordinate for point P is interpolated using the following equation: P v = 
PR(P lv ) + (l-PR)(P 2v ). 

Fig. 3 is a block diagram of an interpolation system that includes an interpolation 
engine 300 in accordance with an embodiment of the present invention. A plurality of 
agents A„ A 2 ,... A„ receives signals 252, 254 and 256, respectively. The signals 
received by the agents are signals from the API and may include information about each 
vertices of a primitive and independent variable X (or Y) from the API.. In one 
embodiment, the agents A 1? A 2 , - A^ may be implemented by pixel calculation unit 
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fa ^362and-span setup unit 366, as shown in Figure 4. The pixel calculation unit 362 

receives, for example, the X coordinates P 1X and P 2X of points Pi and P 2 , respectively, on 
primitive edges (Fig. 2D). For each P 1X and P 2X pair, the pixel calculation unit 362 
generates a request to interpolation engine 300 to calculate a P x value for a point P 
5 within the primitive edges. The span setup unit 366 receives vertex channel information 
for a primitive and a value of an independent variable such as X or Y. The API provides 
the independent variable value. For each independent variable, the span setup unit 366 
generates a request to interpolation engine 300 to calculate X coordinates P 1X and P 2X of 
O points on primitive edges. 

h. f 1 

^ w ft It is further noted that the agents ^Vi, hi, Ah are scalable. An arbiter 258 

%5 receives the outputs of agents Ai, A2,. . . A n and implements a conventional 

fU arbitration algorithm. The arbiter 258 requests the interpolation engine to 

Q interpolate color or texture coordinates U,V and calculate ratio values (e.g., any 

|I of ratio values A' to F' above). The output of the arbiter 258 is received by the 

3 15 interpolation engine 300 which rasterizes primitives by use of direct 

interpolation in accordance with the present invention. The output of router 

260 may feedback into the inputs of agents Ai, A2, A n for purposes of rasterizing 

other points as needed in the primitive. In particular, the ratio values calculated 

Ao A* Art 

f\ by interpolation engine 300 may be routed by router 260 to the agents - Ai, A±, Ah 
20 for use during the rasterization other points in the yrimitivc.Fig rS is a block 
diagram of an interpolation engine 300 in accordance with an embodiment of 
the present invention. The interpolation engine 300 includes a ratio generator 
262, a blending unit 264, and a Z-pipe 266. The output of the ratio generator 262 
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serves as an input signal for the blending unit 264. The blending unit 264 
outputs interpolated values for U and V, while the Z-pipe 266 outputs 
interpolated values for Z. The details and functionality of the ratio generator 
262, blending unit 264, and Z-pipe 266 are discussed further below. 

Referring now to Fig. 6, there is shown an example of one possible 
implementation of an interpolation engine 300 for use in the present invention. 
Field lengths and other descriptive elements are exemplary only, and other 
values may be used in lieu of those shown. From the inputs described above, 
and using the logic described in the above code, interpolation engine 300 
generates interpolated values for U, V, and Ze. Ze represents a screen-based Z- 
coordinate, which may be equivalent to an inverted eye-z value. The ratio pipe 
or engine 262 generates a ratio value R which is also received by the blending 
pipe 264. 

Values for U, V, R and Ze are determined as follows. In one 
embodiment, the ratio pipe 262 includes subtracters 301 and 302 that are used to 
determine the differences between Xi and X, and Xo and X, respectively. As 

stated above, X is an independent variable value such as the physical position at 
the point to be processed. Perspective correction is applied to differences by 
multiplying the output of subtracter 302 by Zo and the output of subtracter 301 
by Zi, using multipliers 305 and 306 respectively. Zo * s the Z-coordinate at Xo, 
while Zi is the Z-coordinate at Xi. 

An adder 307 adds the outputs of multipliers 305 and 306, and a module 
308 inverts the output of adder 307 to generate the term 1/ [ z 0 (X- Xo) + Zi(Xi- 

X)]. A multiplier 309 then multiplies the outputs of module 308 and multiplier 
306 to generate the ratio R = Zi(Xi-X)/[ Zo(X- X 0 ) + Zi(Xi-X)]. 
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In one embodiment, the Z pipe 266 includes a subtractor 315 that 
determines the difference between Xi and Xo. A multiplier 317 multiplies the 
output of subtracters 315 and module 308 to generate Ze = (Xi- Xo)/ [Zo(X- Xo) + 
Zi(Xi-X)]. The Ze value is a perspective depth coordinate that represents the 
inverse value 1/w wherein w is the above mentioned perspective coordinate 
parameter. It is noted that the Z pipe 266 is an optional element in the block 
diagram of Fig. 6. 

In one embodiment, the blending pipe 264 includes subtracters 320 and 
322 that are used to determine the differences between Uo and Ui, and Vi and 
Vo, respectively. A multiplier 324 multiples the outputs of subtractor 320 and 
multiplier 309 to generate the output R(Uo- Ui). An adder 326 adds the output 
of multiplier 324 with Ui to generate U = R(Uo- Ui) + Ui = R Uo + (1-R) Ui. 

A multiplier 328 multiplies the output of subtractor 322 with R to 
generate R(Vo- Vi). An adder 330 adds the output of multiplier 328 with Vi to 
generate V = Vi +R(Vo- Vi) = R Vo + (1-R) Vi. 

The elements in the blending pipe 264 may be scalable, as represented by 
dot symbols 335. In addition, based on the teachings of the present invention, 
one of ordinary skill in the art will realize the mathematical functions provided 
by the blending pipe 264 may be varied. As a result, the elements of blending 
pipe 264 (e.g., multiplier 324 and 328 and adders 326 and 330) may be connected 
in other configurations to achieve other desired mathematical functions for the 
blending pipe 264. 

The interpolated output values U,V for perspective interpolation are 
listed and described further below. 
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U = Ui + [Zi(Xi - X)] (Uo -Ui)/ [(Xi - X) Zi + (X -Xo) Zo] 
(Eq.4) 

V = Vi + [Zi(Xi - X) ] (Vo -Vi)/ [(Xi - X) Zi + (X-Xo) Zo] (Eq. 5) 

The term (X -Xo) Zo * s °* a small value and is a term that prevents, for 
example, the ratio Zo(Xi - X)/[(Xi - X) Zi + (X-X 0 ) Zo] from having a unity value. 
The perspective coordinate parameters Zi and Zo are applied to the terms in 
equations 4 and 5 above to the interpolated values for U,V. 

As stated above, the blending pipe 264 if Fig. 6 is scalable. Hence 
subtracter stages may be added to calculate, for example differences in X 
coordinates (Xo - Xi) and/ or differences in Y coordinates (Yo - Yi) . A J ^rr^b % 
multiplexer stage may be implemented in the interpolatioryofTig. 6 so that 
inputs can be selected among the values such as Uo, Ui, Vo, Vi, Xo, Xi, Yo, Yi or 
differences between these values (e.g., Xo - Xi). It is also noted that if multiplier 
stage 305 and 306 are omitted, then a linear interpolation function may be 
performed by the interpolation engine. Alternatively, a linear interpolation 
function may be performed by suitably blocking the Zo and Zi inputs. 

Subtracters 301, 302, 303, 310, 319, 320, and 321, and adders 333, 334, 335 
can be implemented as conventional floating-point, integer or other subtracters, 
as specified. Subtracters are known by those of ordinary skill in the art, and the 
type of subtracter chosen is not intended to be limiting. Alternatively, one 
embodiment of the present invention employs an improved subtracter as 
disclosed in related application Serial No. 09/020,152 for "Fast 
Adder/ Subtracter for Signed Floating Point Numbers", by Konstantine Iourcha, 
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Andrea Nguyen, and Daniel Hung, filed February 6, 1998, the disclosure of 
which is incorporated herein by reference. 

It is within the scope of the present invention to implement a program or 
code that can be stored in an electronically-readable medium to permit a 
5 computer to perform the methods described above. 

From the above description, it will be apparent that the invention dis- 
closed herein provides a novel and advantageous system and method of graphic 
interpolation for rendering primitives. The foregoing discussion discloses and 
describes merely exemplary methods and embodiments of the present 
n w invention. As will be understood by those familiar with the art, the invention 
pTj may be embodied in other specific forms without departing from the spirit or 

1 1 essential characteristics thereof. In particular, other component types may be 

^ employed, and other primitive types, dimensions, or projections may be used. 

= y Accordingly, the disclosure of the present invention and example figures are 

5 ^ intended to be illustrative, but not limiting, of the scope of the invention, which 

|y 

H is set forth in the following claims. 
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